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I. INTRODUCTION 

Field theories in two spatial dimensions have long been 
recognized as important for understanding several phys- 
ical phenomena that can be well approximated as planar 
ones, such as high temperature superconductors and the 
fractional quantum Hall effect. In the last years there 
has been an enormous interest in these two dimensional 
systems in condensed matter physics. In particular, a 
large number of important new phenomena discovered 
recently in condensed matter lies in this class. As exam- 
ples we can cite the metal-insulator transition (MIT) [1| 
and graphene, an isolated single atomic layer of carbon, 
in which electron transport is essentially governed by 
Dirac's relativistic equation 0. It is then feasible to 
treat many of these planar (or quasi planar) systems in 
condensed matter in terms of quantum field theory mod- 
els for fermions in 2-1-1 dimensions (two space like and 
one time like coordinates). 

Among the many field theory models useful for un- 
derstanding a plethora of phenomena in condensed mat- 
ter systems, those that include a four-fermion interac- 
tion have been extensively used in the context of stud- 
ies of planar systems. In particular, Nambu-Jona-Lasino 
type of models Q in 2-1-1 dimensions and, among them, 
the Gross-Neveu (GN) model jj], have been employed 
to study these systems. For example, the GN model has 
been considered as an appropriate model to study low en- 
ergy excitations of high temperature superconductors ^5J , 
while analogous models with generic quartic fermionic in- 
teractions have also been used to study quantum prop- 
erties of graphene @ . 

In these planar systems, GN field theory models 
for fermionic interactions are commonly used to study 
their symmetry properties, in particular chiral symme- 
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try breaking and restoration, at finite temperature and 
densities and also in the presence of external (magnetic) 
fields. In this paper we will be interested in investigat- 
ing the effects of how an "asymmetrical doping" can af- 
fect the chiral symmetry, or the metal-insulator transi- 
tion in a GN type of model for a two-dimensional sys- 
tem of fermions. An asymmetrical doping can be defined 
as an imbalance between the chemical potentials of the 
electrons with the two possible spin orientations ("up" 
=1, and "down" =|) inserted in the system by a dop- 
ing process. Since the densities of the t and J, electrons 
are directly proportional to their chemical potentials, an 
asymmetrical doping is equivalent to an asymmetry in 
the chemical potentials for spin up and spin down elec- 
trons. This chemical potential asymmetry can be pro- 
duced by the effect of an external magnetic field parallel 
to the system's two-dimensional plane. 

A magnetic field applied parallel to the system's two- 
dimensional plane couples only to the spin of the elec- 
trons, but not to the electrons orbital motion. Therefore, 
Landau levels that would appear due to the coupling of 
a (perpendicular to the plane) magnetic field to the elec- 
trons' orbital motion do not appear here. Instead, the 
in-plane magnetic field generates an intrinsic Zeeman ef- 
fect which polarizes the system. This is because at zero 
magnetic field electrons with spin up and spin down have 
the same density, but an in-plane magnetic field, due to 
the Zeeman effect, causes a difference between the spin 
up and spin down densities. Recent studies in the con- 
text of graphene 0, Q have suggested that the Zeeman 
effect can be very important for the electronic proper- 
ties of these systems. The intrinsic Zeeman effect is thus 
useful to reveal the important role played by the spin 
degree of freedom of the electron and the polarization of 
the system 0. 

Here we will focus on the properties of a planar system 
modeled by a GN four- Fermi interacting model and study 
how a spin density asymmetry influences on the chiral 
symmetry of the system, i.e. on the role of the Zeeman 
contribution to the system's magnetization. The rest of 
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this paper is organized as follows. In Sec. II, we review 
the GN model and the effect of including an in-plane 
constant magnetic field and the Zeeman effect. In Sec. 
Ill we evaluate the Landau's free energy (the effective 
potential) of the model and determine how the chiral 
symmetry is affected by the spin density asymmetry. In 
Sec. IV we determine the magnetic properties of the 
system. Finally, in Sec. V we present our concluding 
remarks. 



II. THE MODEL ACTION 



For example, in one space dimension, Eq. ([2]) is the 
equivalent of the Takayama-Lin-Liu-Maki (TLM) model 
the continuous version of the model proposed by Su, 
Shrieffer, and Heeger (SSH) for polyacetylene [H, [Tj], in 
the adiabatic approximation, where lattice vibrations are 
neglected, and used to study the metal-insulator transi- 
tions in general (see e.g. [l4| and references therein). 

Let us now consider the application of a generic ex- 
ternal magnetic field to the system and its effects. It is 
convenient that to start by writing the grand canonical 
partition function for the Lagrangian density model 



We start by considering a planar four- Fermi model de- 
scribing interacting fermions with a Lagrangian density 
given by 



s=T,i 



(1) 



where tp'^ is a four-fermion field with N flavors and s is 
an internal symmetry index (spin) that determines the 
effective degeneracy of the fermions. In Eq. ^ a sum 
in the flavors is implicit. For example, for materials such 
as polyacetylene or graphene, N = 2. We will keep TV 
general throughout this work for convenience. 

In Eq. ID), A is a coupling term and wf is the Fermi ve- 
locity. The gamma matrices are 4x4 matrices and we fol- 
low the representation given e.g. in Ref. lOj for fermions 
in 2-1-1 dimensions. In this case, the model H]) possesses 
a discrete chiral symmetry, ijj — s- 75'0, ip — > — i/'Ts? with 
the 75 matrix defined as in [lOj . This discrete chiral sym- 
metry is the one considered in this paper, along with its 
breaking and restoration conditions. Note that it is bro- 
ken when a gap (or a nonvanishing vacuum expectation 
value for (iptp)) is generated. 

It is useful to rewrite the fermion interaction term in 
Eq. Il]) in terms of a boson field A, in which case ii[4',ip] 
becomes 



N 



ihdt — ihvFj-^ — A 



r 



(2) 



A can be seen, for example, as representing the coupling 
of electrons to a local value of the dimerization. Equa- 
tion 12]) is equivalent to Eq. H]) as can be easily verified 
by using the field equation for A(a;) and substituting it 
back into the Lagrangian density, re-obtaining the char- 
acteristic four-Fermi interaction. 

Models of the type of Eq. ^ , or equivalently Eq. ^ , 
are of the form of a Gross-Neveu model [3| , with applica- 
tions found in many areas of condensed matter physics. 



Z = J DAY[d,P^D^ ei^v{-SE[^,^,A]} , (3) 



where the Euclidean action SE['4'y'4']j from the La- 
grangian density Eq. reads 



h0 



A] = / dr I d'x {Y, r [ndr 



-I- A + 70 + Ylogf^BB\\ 



where /? = l/{kBT), ks is the Boltzmann constant, /i is 
the chemical potential. Ax and Ay are the vector poten- 
tial components ( e.g. corresponding to a magnetic field 
perpendicular do the system's plane), B\\ is the magnetic 
field parallel to the system's plane and crs7o .g /i_BS||/2 
is the corresponding Zeeman energy term, with (j-^ — 1, 
CTj = — 1, 5 is the spectroscopic Lande factor and //s is 
the Bohr magneton. We must point out that graphene 
samples have been recently studied in strong magnetic 
fields (up to 45 T) In these experiments the mea- 

sured effective g* factors were found very close to the 
bare electron g factor {g = 2). In our work developed 
below, we will not assign a specific value for g and we 
keep it also general (like N) for convenience. 

By choosing a gauge where the three-dimensional vec- 
tor potential is given, for example, hy A = (0, B±x, B^), 
we see from Eq. (U) that B± couples to the orbital mo- 
tion of the electrons and it will result in the Landau levels 
for the system in this magnetic field. The parallel mag- 
netic field couples to the electrons' spin and produces the 
Zeeman energy term in Eq. ([?]). From the form of the 
Zeeman energy term in Eq. ([4]) we see that it can be 
added to the chemical potential, thus defining an effec- 
tive chemical potential term in the action of the form. 



(5) 
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where /i| = ^+S^, and = n—Sfj,, with 5^ = g ^bB\\/2. 
The role of /i can be interpreted as to account for the ex- 
tra density of electrons that is supplied to the system by 
the dopants, while (5/i measures the amount of asymme- 
try introduced and it is directly proportional to the in- 
plane applied external magnetic field. As we explained in 
the introduction, in this work we will be concerned with 
the effects of the Zeeman term, so from now on we take 
= and only consider a constant in-plane external 
magnetic field Bv. = Bp- 



for Ac for studying how an asymmetry between the up 
and down fermions' spins, generated by the constant in- 
plane magnetic field, will change the phase diagram of the 
model. Possible phenomenological applications of these 
results for systems like for example graphene and pla- 
nar organic conductors, as polyacetylene, will then be 
discussed. 

The effective potential is defined from the grand canon- 
ical partition function ([3]) by 



III. THE SYSTEM'S EFFECTIVE POTENTIAL 
UNDER EXTERNAL EFFECTS 

In the applications with models of the GN type, we 
are mostly interested in studying the effective potential, 
or Landau's free energy, for a constant scalar field con- 
figuration Ac, in which case chiral symmetry breaking 
and dynamical fermion mass generation can be most con- 
veniently studied. Here, we use the effective potential 



Foff = -— InZ 
/3V 



(6) 



where V is the volume. Then, from the partition func- 
tion Eq. ([3]) with (dl and using Eq. ([5]), for a constant 
background scalar field Ac and in the mean-field approx- 
imation, which is the same as considering just the leading 
terms in a 1/N expansion, or the large- approximation 
obtain that the effective potential VcS is given 

by 



J 



v;ft(Ac,r,/it,/i|) 



N 



2hvFX 



+ 00 „ 

Al-NksT J2tr Y: J 



s=T.i n=-oo' 



r 



j^--^\n[{~iujn + ^is) -VFjoJ.p~-foAc] , (7) 



where ujn = (2n+ Vj'KkBT are the Matsubara frequencies 
for fermions. Note that in the absence of asymmetrical 
doping, = ^1 = fi, or 6^ = 0, Eq. ([7]) just repro- 
duces (in the natural units where h = — vp ~ 1) 
the effective potential of the Gross-Neveu model in 2+1 
dimensions 

Performing the sum over the Matsubara frequencies 
and taking the trace in Eq. ([7]), we find 



Ve«{Ac,T,fl^,fll) 



N 



2hvpX 



-A 



-iln(l-|-e-'^^f)+iln(l + e-^^r 



■ In 1 
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where E^^ — Ep±fi^^i and Ep = y/ v^pp"^ + A^. In terms 
of the band structure, we can interpret the result seen 
in Eq. ([5]) as like 6^ has lifted the degeneracy of the 
conduction and valence bands in the matter part of V^s- 
At zero temperature and chemical potential and in the 
absence of the external magnetic field, Kff becomes 



Feff(Ac) 



A^ 



2UvfX 



-At - 2N 



d^p 

{2TTny 



+ A2 . (9) 



By using a momentum cutoff A to regulate the vacuum 
divergent term of Vcs{Ac) and by defining a renormalized 
coupling Xfj as 



J_ ^ hvF d^V^sjA,) 



N 



dAl 



(10) 



where Aq is a renormalization point, that can be cho- 
sen, as usual, as given by the nontrivial minimum of the 
renormalized effective potential. 

In terms of Xh, the renormalized effective potential 
reads (after subtracting an irrelevant field independent 
divergent vacuum term) 



Veff,R(Ac) 



N 



2hvFXR 
N 

'K{hVFY 



I Ac 



|Ao|A; 



(11) 



The nontrivial minimum of Vef f.R{Ac) can now be easily 
found and it is given by 



Ao 



Tivft^ 



(12) 



At Ac = Aq, the effective potential reads 
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Voff,R(Ac = Ao) 



N 



A3 



(13) 



which shows that Kff(Ac Ao) < T4//(Ac = 0) = 0. 
The non-trivial solution is then energetically preferable 
for the (undoped) system, which then corresponds to a 
(dynamical) gap, i.e., the presence of a chiral nonvanish- 
ing vacuum expectation value. 

The effective potential at finite chemical potentials and 
in the zero temperature limit, from Eq. ([5]), is given by 



V;ff,R(Ac,Ait,x) = Kff,R(Ac) 

T 

Jo Znti 

r t^(E,-\f,,\) , (14) 

where V;ff,R(A,) is given by Eq. Bi^ = e{fil^-Al) 
is a step function, defined as Q{x) = 0, for cc < 0, and 
Q{x) = 1, for X > 0, and is the Fermi momentum of 
the f (J.) fermions, 



1 



(15) 



By performing the momentum integration in Eq. (|14p . 
we obtain 



Voff,R(Ac,Ait,x) = T4ff,R(Ac 



47r(?iuF)^ 
ATr{hvFy 




Minimizing T4ff,R(Ac, /X|_|) with respect to Ac yields 
again the trivial solution (|Ac| = 0) and 



|A,|-Ao + ^(a/T-|Ac|) 



(K|-|Ac|) = 



(17) 



Before continuing, lets us specialize to the symmetric 
limit, (5/i = 0. In this case Eq. pT]) reads |Ac| — Ao + 
8(/x^ - A^) (^ - |Ac|) = 0. To solve this equation we 
need to know the critical chemical potential fic at which 
the symmetry is restored. This is found through the 
equation Vcs^b.{Ac = 0, ^lc) = T4fi,R(Ac = Ao, /Ltc), which 
yields fic = Aq. Thus, the ground state of the symmetric 
system is characterized by 



At and beyond the doping level fj,c, we note that the 
energy of the fermions is linear, Ep = vpp, characteristic 
of two dimensional gapless systems [isj . 

Considering now the asymmetrical system case, we see 
that by applying a constant magnetic field Bq parallel 
to the plane, the asymmetry 5 pi is increased as in- 
creases. We study the asymmetry relative to ^, and look 
at the system's ground state for a given asymmetry. In 
Fig. [T] we show the effective potential Kff,R(A, /i|,|) as a 
function of A^ for various asymmetries < dfj,/ Aq < 1. 
The first curve in Fig. [T] is the effective potential Vcs^r 
as a function of Ac/Aq for ^ — O.Sfic and Sfi/Ao — 0, 
with the minimum of Kjff^ at Ac = Ao. The second 
curve from top to bottom shows the symmetry restora- 
tion, where fj, = fj,c, and Sii/Aq = 0. The following curves 
are for Sji/ Aq = 0.5 and (5/i/Ao = 1.0. These results, in 
particular, show that there is not a value for an asym- 
metry (5/i for which a new minimum, where Ac{5^) ^ 0, 
could appear. This is clearly seen from the third and 
fourth curves in Fig. [T] (from top to bottom), which show 
that the minimum of the effective potential is always at 
Ac = 0. This result should be contrasted with the one 
obtained for the one-dimensional space system case re- 
cently studied p^ . where by increasing the asymmetry, 
a new minimum for the effective potential emerges at a 
critical chemical potential asymmetry 5^c- Thus, a non- 
vanishing and stable gap beyond (5/ic would exist, till it 
disappears completely at 3> (J/ic, through a second or- 
der phase transition. The absence of this new minimum 
in our case will be explained below. 

Finally, let us now consider the complete renormalized 
effective potential at both finite chemical potential and 
temperature, obtained from Eq. ([5]) after performing the 
momentum integrals. It is given by 



yoff(Ac,T,^l,^l) = T4ff,R(Ac) 



Acl . 



,-/3(Ae-A'T)l 



+ Li2[-e-'^(^=+''T)]} 
- ^{Li3[-e-'^(--T) 



Lis 



} + (19) 



where lAu{z) is the polylogarithm function and it is de- 
fined (for J/ > 0) as [20| 

oo ^ 

fc=i 

From Eq. (jl9p we can now verify how the Zeeman 
term, manifested by the density asymmetry term dji, 
changes the usual (for djj, = 0) chiral phase transition in 
the GN model when both finite chemical potential and 
temperature are considered. We start by deriving the 
gap equation. 



Ac = 



Ao, for < pLc 
0, for fi > fJ-c 



(18) 



d 
dA, 



V;ft(Ac,T,/i|,^x) 



(20) 
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FIG. 1: The effective potential at zero temperature, Eq. (|16p . 
in units of NAQ/{hvF)^- The top curve (solid line) is for 
fj, = 0.8/ic and (5/i/Ao — 0, with the minimum of Vcff.R at 
Ac — Aq. The dotted line is for /i — fic, and Sh/Aq = 0, 
with the minimum of V^ff^R at Ac — Aq and at Ac = 0. 
The following curves are for Sjj,/ Aq — 0.5 (dashed line) and 
5/i/Ao — 1.0 (dash-dotted line). 
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FIG. 2: The critical curve Ac(T, fi, Sjj,) — for different values 
of asymmetry. The top curve (solid line) is for Sfi = 0, while 
the other curves are for 5fi — O.S/^c (dotted line), Sjj. — 0.8/ic 
(dashed line) and 5fi = 0.95^c (dash-dotted line), respec- 
tively. 



" ■ * % , 



^ ^ 




- — - 


\ '• \ 

\ \ 

\ -A 




\ •.\ 










\ •. 












X \ 

















which gives 

Ac = Ao- 
+ hi 

+ hi 



1 



|ln 1 + p-l^i^c+f^+Sfj^) 



2/3 

I _|_ g-/5(Ac-A'-t)7') 

1 + e"'^('^<="l''"*'''l) 



In 



^-l3{A^+\fj.-5fj.[ 



(21) 



It can be easily checked that the T = hmit of the above 
equation reproduces the result (|17p . 

The critical curve Ac{T, ^,dfi) — 0, that is obtained 
from Eq. ([2T|) in the case 6fi = 0, just reproduces the 
known result [lO|, with a line for a second order phase 
transition in the fi — T plane, starting at the critical point 
= 0,T = Tc), where ksTc = Ao/(21n2), and ending 
in a first order critical point at (/i = fJ-c,T = 0). By 
increasing the asymmetry (the magnitude of the parallel 
magnetic field Bq) the effect is to promote chiral symme- 
try restoration, as can be explicitly seen in Fig. [2l 

We can now interpret the nonexistence here of new 
minima at /i = /Xc in the presence of an asymmetry, in 
contrast to the findings of Ref. 19] for the one space di- 
mension case. This can be traced to the nonexistence 
of a critical line for first order chiral phase transition in 
the GN model in two space dimensions in the mean-field 
approximation. The phase diagram of the GN model in 



one space dimension, in the mean-field approximation, 
has a second order critical transition line in the jJ, — T 
plane that meets a first order transition line at a tricrit- 
ical point. This is actually a typical phase diagram seen 
in many other four-Fermi interacting models, including 
Nambu- Jona-Lasino type of models in three space dimen- 
sions. However, this typical phase diagram was absent in 
the two space dimensions GN model until recently [2l| . 
In Ref. [2l| it was shown that terms contributing to the 
effective potential going beyond the mean-field approxi- 
mation would produce a first order critical line merging 
with the second order critical line at a tricritical point. 
The effect of the asymmetry close to the critical point 
/ij, seen in Ref. [l^ would then probe the metastable re- 
gion around the first order critical line, manifested by the 
formation of a new local minimum in the presence of an 
asymmetry. The same should be seen here, if corrections 
beyond the mean-field approximation would have been 
considered. We do not see this formation of a new local 
minimum close to the first order critical point here be- 
cause the mean-field approximation in the two-space di- 
mensional GN model misses this metastable region. We 
expect, though, that these corrections would produce a 
very small effect here. This is because the metastable re- 
gion seen in Ref. [2l[ was considerable smaller than the 
one in the one-space dimensional GN model, so its con- 
tributions to our derivations should, likewise, be small. 



6 



IV. MAGNETIC PROPERTIES 



As we have already seen, the imbalance in the chemi- 
cal potentials of the t and J, electrons are induced by the 
application of a static in-plane magnetic field in the sys- 
tem, with a Zeeman energy Ez = ±gfJ-BBo/2 and, in the 
present case, 6fi — \Ez\ — gfiBBQ/2. As a consequence, 
the number densities and n^, defined by 



(22) 



are obviously imbalanced due to the asymmetry between 
/i| and /ij^, and will depend on dfi. Likewise, the spin 
polarization of the system as a result of the Zeeman ef- 
fect produces a net (Pauli) magnetization of the system, 
which is defined by [24I 



n 



n 




and a magnetic (Pauli) susceptibility. 



(23) 



(24) 



Note that the polarization state of the system is deter- 
mined by the asymmetry, which depends on the intensity 
of the applied magnetic field. As we change the magnetic 
field, and then the polarization, we expect the magnetic 
properties of the system, e.g. the magnetization and the 
susceptibility, to change as well. As we are going to see, 
the change in behavior can be quite remarkable, mostly 
when the susceptibility as a function of temperature is 
concerned. 

Let us start by presenting some results for the number 
densities for fermions of spin up and spin down, as given 
by Eq. (|22p . We consider first the case oi fj, = which 
represents the undoped system, i.e., the insulating state 
at which Ac = Aq. Since /i = we have effective chem- 
ical potentials given by /i| = Sfi and = —6fi. From 
Eq. ((22| and the results of the previous section, one sees 
that if Sfi < Aq, than = rt| = 0, giving M — x = 0. 
On the other hand, if Sfi > Aq, there are non- vanishing 
densities, but they are always equal (n^ = ni), resulting 
again in M = x = 0. This means that the undoped (in- 
sulating) system, for which /i = 0, is never magnetized 
(polarized) at any temperature. So, let us then consider 
the cases where /i 7^ 0. In Fig. [3] we show how the spin 
up and down fermion densities (in units of iVAg/(?iui?)^) 
change with the asymmetry (i.e., when the applied mag- 
netic field increases), for a chemical potential chosen here, 
as an example, as fi — 0.7/ic = O.TAq and for tempera- 
tures ksT — O.5A0 (solid and dashed lines, for and n^, 
respectively) and for ksT = Aq (dash-dotted and dotted 
lines, for n| and n^, respectively). Note that in both 
cases ni initially decreases with the increasing magnetic 







FIG. 3: The number densities of spin up and down fermions 
as a function of the asymmetry, for /i = 0.7Ao and for tem- 
peratures ksT — O.5A0 (solid and dashed lines, for and 
rij, respectively) and for ksT = Aq (dash-dotted and dotted 
lines, for n| and n^, respectively). 



field, vanishing at the full polarization point /i = 6fi, as 
expected, and starts increasing again beyond that point. 
This is when the magnetic field becomes strong enough 
to be energetically favorable to spin down fermions of 
the filled valence band to be promoted to the conduction 
band. This occurs only at finite temperature, since at 
zero temperature Udown is always zero at and beyond 
the critical magnetic field Be (or the full polarization 
point) [1^. Note that the spin up fermion density al- 
ways increase with the applied magnetic field (which is 
in the same direction of the spin up fermions). The plots 
also indicate that the rate of spin down fermions turning 
into spin up ones (below the full polarization point) as 
the result of increasing the applied magnetic field is con- 
stant and is approximately the same rate of promotion 
(above the full polarization point) of spin down fermions 
of the valence band to the conduction band. Another re- 
sult we can notice from Fig.[3]is that the difference of the 
two densities, i.e. the Pauli magnetization Eq. (|23p . is 
always positive with the increasing magnetic field. This 
is true for any other values of chemical potential and 
temperature. 

Next, in Fig. 2] we present the results for 
the Pauli magnetization, expressed in units of 
N {guB/ 2) Aq/ (Tivp)'^ , for different values of tem- 
perature (and again fixing the chemical potential as 
fi = O.7A0 for comparison purposes) when the asym- 
metry is increased (or, equivalently, in terms of the 
Zeeman field Bq). Here we can easily see a clear change 
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FIG. 4: The Pauli magnetization as a function of the asymme- 
try, for fi = 0.7Ao and for ksT = 0.5Ao (dotted hne), 0.7Ao 
(dash-dotted hne), 0.9Ao (dashed line) and l.lAo (sohd hne). 



in behavior of the magnetization as the full polarization 
point is crossed. The most important observation we can 
notice from the results seen in the Fig. [4l it is that below 
the full polarization point, here given by dfi = 0.7/ic 
(which is the value of /i that we have considered in 
the analysis), the larger is the temperature, the larger 
is the rate of increase in the magnetization, while has 
an opposite behavior above the full polarization point. 
This then will reflect remarkably in the Pauli magnetic 
susceptibility results. As the temperature is increased, 
the Pauli magnetic susceptibility should as well increase 
for magnetic fields below of that which gives the full 
polarization. Above the full polarization magnetic field, 
the magnetic susceptibility should decrease with the 
increasing temperature. This is confirmed by the results 
presented in Figs. [H] and O for the cases below the full 
polarization point and at and above it, respectively. 

Figures [5] and El show the magnetic susceptibility, ex- 
pressed in units of N{gfiB Aq/ (hvp)'^ , as a func- 
tion of temperature, when the asymmetry is varied from 
dfi ~ O.lAo up to l.lAo, with fixed chemical potential 
/i = 0.7Ao. For all other values of chemical potential we 
find similar behavior for x{T, fJ-, Sfi). 

Analyzing the behavior of the magnetic susceptibility 
X as a function of the temperature, we find that, for 
those cases below the full polarization point, e.g. Fig. O 
in the chiral broken (insulating) phase (Ac ^ 0), x is 
nonlinear for T < Tc, where the values of Tc are indicated 
by the horizontal dotted line in Fig. [51 This horizontal 
critical line can be determined explicitly and it is given 
by N{gpLBY Ao/[47r(?iui?)^]. When this critical line is 



FIG. 5: The Pauh magnetic susceptibility as a function of the 
temperature for /i = 0.7Ao, for the values of the asymmetry 
5^l (from bottom to top) 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 (in units 
of Ao). The dotted line indicates where the phase transition 
happens. 



crossed and we go from the broken (insulating) phase to 
the symmetric (metal) phase, for T > (or Aj, = 0), 
the magnetic susceptibility becomes a linear function of 
the temperature. Note that the magnetic susceptibility, 
below the full polarization point, is always an increasing 
function with the temperature, as already anticipated by 
the behavior of the magnetization seen from Fig. 2) 

As we go above the full polarization point, the mag- 
netic susceptibility changes again its behavior with the 
temperature. We checked that at high temperatures, 
ksT 3> Aq, it goes exactly as predicted by the Curie- 



Weiss law [221, x(m < (5^,^ > Ao) - 1/T. The turning 



point of behaviors for the magnetic susceptibility, as seen 
in Figs. [5|and[6l is the full polarization point ^ — Sfx. 

A few analytical results for the magnetic susceptibility 
can be obtained explicitly. At zero temperature, for any 
fi + S^J, < /ic, we find that x{T = 0, fi, (5/i) = 0, while for 
fi + S ^j. > fj.c vfe find that xiT = 0, /i, (5/i) 7^ 0. The two 
curves from bottom to top seen in Fig. [5| correspond to 
the former case, where fj, + Sfi < fic- For the cases where 
fi + 5^ = fic, with /i < /Zc or (5/i < /ic, which in Fig. \E\ 
corresponds to the third curve from bottom to top, we 
find the analytical result for the magnetic susceptibility 
at T = 0, 



X(T = 0,fi + S^i = He) 



167r(?iUf^)^ 



Ar 



At /i = /ic, for any (5/i, or fi + Sfi > fic such that Ac 



(25) 
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FIG. 6: The magnetic susceptibility as a function of the tem- 
perature for fi — 0.7Ao, for the values of the asymmetry 
Sfi = 0.7, 0.8, 0.9, 1.0, 1.1 (in units of Aq). 



to the value given by Eq. (|25p . for ^ + 5^ = ^c, with 
fi < and Sfi < /ic, and then it jumps again from the 
value given by Eq. (^5]) to the limiting lower value of sus- 
ceptibility given by Eq. (|27p . obtained for fi = 0.5/ic- For 
any other possible value that the magnetic susceptibility 
can have, it will be larger than the value obtained from 
Eq. (P7|) . This behavior for the magnetic susceptibility 
could in principle be tested experimentally, induced by 
either a change in chemical potential (e.g. by increasing 
the doping concentration in planar systems) or by an in- 
crease in the magnetic Zeeman field (thus increasing the 
spin asymmetry). 

We also note that, as we have seen in the previ- 
ous section, since for any value of /i < /ic, the value 
of asymmetry Sfi = will lead to the phase transi- 
tion. This value of asymmetry corresponds to a criti- 
cal magnetic field for the chiral phase transition. Since 
= Aq — g^BBo_c/'2, we find that this critical mag- 
netic field is given by 



where we have used that /ic = Aq . We may compare this 
result, for example, with that of MIT [H. [23l|: 



(see Fig. [Ij, i.e., in the metal (chiral restored) phase, we 
find that 



x(T = 0,Ac = 0) 



Ao. 



(26) 



This, in particular, corresponds to the values of parame- 
ters determining the critical line seen in Fig. [5l 

The full polarization point can be reached for various 
dopings and in-plane magnetic fields, i.e., always that 
fi = Sfj,, in which case /i| = 0. In particular, at the full 
polarization point, for any /i > 0.5/ic, we find that 



X(T = 0, /i > 0.5/ic, /i = (5/i) 



(27) 



We have also verified that this same result for the sus- 
ceptibility is also obtained for other cases outside the full 
polarization point, e.g. for /f > Mo with any value of the 
asymmetry Sfi, or for /i > 0.5/ic and Sfj, > fic- Note that 
this result includes the case leading to Eq. for the 

special case of /i = /ic. The curves starting at the dot- 
ted line in Fig. [5] correspond exactly to examples of these 
cases, with the dotted line obtained when /i = 0.7Ao is 
substituted in Eq. (gT]). 

The parameters determining the result given by 
Eq. ([TT]) , together with those determining Eq. (^5]) , show 
that there are two major regions of parameters, at zero 
temperature, that preclude the Pauli magnetic suscep- 
tibility for attaining a value. We find that the Pauli 
magnetic susceptibility jumps discontinuously from zero 



Bpoi — 



2Ef 



(29) 



where Ep is the Fermi energy, g* is the effective Lande 
g factor and m* is the effective mass. 



V. CONCLUSIONS 

To summarize, we have investigated the mean- field 
phase diagram of planar systems, that can be modeled 
with a four-Fermi type of model, upon asymmetric dop- 
ing. We have obtained the magnetization and magnetic 
susceptibility for this system. The analysis made in this 
work can be of relevance in studies of many planar sys- 
tems of interest in condensed matter, such as, for ex- 
ample, organic conductors made of polymer films and 
graphene. In particular, regarding the zero temperature 
magnetic susceptibility of these systems, we have pre- 
dicted that it can change abruptly from zero to the value 
given by Eq. (j25[) . when increasing either the doping or 
the magnetic Zeeman field, such that fj,^ crosses a criti- 
cal doping value given by /ic. Then, it can jump again 
abruptly up to the point given by Eq. (j27p . with the 
minimum value attained for /i — 0.5/ic. This is a direct 
prediction of our results that could be tested experimen- 
tally 

As far the dependence on temperature is concerned, we 
have obtained that the magnetic susceptibility can have 
quite different behaviors depending whether the system 
is below or above the full polarization point /i = Sfi. Be- 
low the full polarization point, fi < Sfi, while in the chiral 
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broken phase the magnetic susceptibihty increases non- 
hnearly with the temperature, in the chiral symmetric 
phase it tends to a hnear function of the temperature. 
Unfortunately, it seems that there is still so few or none, 
to our knowledge, experimental results for the magnetic 
susceptibility of planar systems of the type that could 
be modeled by the model studied here, such as planar 
polymers or graphene, so to be able to compare with our 
findings. Even so, despite the very different nature, we 
recall that there are measurements of the behavior of the 
magnetic susceptibility with the temperature for metal 
alloys, most notably titanium based that displays 
exactly the same behavior as we found here, increasing 
almost linearly with the temperature, while the sudden 
change in the susceptibility at the chiral phase transition 
point seems to be analogous to the behavior seen in the 
measured susceptibility of blue bronze (K0.3M0O3) [2^ 
and indicative of a Peierls transition there. The behav- 
ior of the magnetic susceptibility changes again as we 
go above the full polarization point. For ^ < Sfj, the 
magnetic susceptibility decreases when the temperature 
increases. In particular, for high temperatures, T ^ Aq, 
it obeys the Curie- Weiss law. The applied magnetic field 
is then determinant on the type of the behavior observed 
for the magnetic susceptibility. There is a critical mag- 
netic field, proportional to the chemical potential (the 
doping concentration) that determines the two behaviors 
for the magnetic susceptibility as a function of the tem- 
perature. We recall that there is a similar behavior seen 
in pristine graphite [2^ . where for low magnetic fields 
the magnetic susceptibility is observed to increase with 
the temperature, while for large fields it goes down with 
the temperature. The major difference there with what 
we see here is that graphite is diamagnetic, while our 
results reflect the behavior of a paramagnetic material. 
We hope that in the future there will be measurements 
of the Pauli magnetic susceptibility for planar systems so 
to be able to more closely compare with the results we 
have obtained here. 

Our results have also shown that the Zeeman effect 
in these planar systems tends to weaken the chiral sym- 
metry, thus, the insulating to metal transition may hap- 
pen at a smaller critical temperature in the presence of a 
Zeeman field. This behavior, due to an increasing mag- 



netic Zeeman field, is exactly the opposite to what is 
observed when a perpendicular magnetic field is applied 
to these systems (see e.g. [13, Ull), in which case the chi- 
ral symmetry breaking becomes stronger by the effect of 
the magnetic field and, thus, the transition happens at 
a higher temperature in the presence of a perpendicular 
magnetic field. These two opposite effects caused by mag- 
netic fields, when applied parallel or perpendicular to the 
system's plane, can be an useful tool to regulate the in- 
sulating/metal behaviors for these planar systems, when 
the parallel and perpendicular fields are applied simul- 
taneously and independently, opening interesting possi- 
bilities for uses of these type of materials in practical 
applications as electronic devices. Further studies on the 
magnetic properties of these systems we hope to make 
and to report on them elsewhere in the future. 

Finally, since in most realistic experiments the mea- 
sured quantities of interest are related to electric trans- 
port, it is expected also in those cases the Zeeman split- 
ting to have important effects. However, a calculation of 
conductivity effects and transport properties cannot be 
addressed with the methods we used here, based on the 
calculation of the effective potential (free energy), which 
are more suitable for the analysis of the phase structure 
of the model. But based on experimental studies of an 
in-planc magnetic field, for example in graphene, the Zee- 
man splitting has been shown to be important in both 
the spin transport and conductance fluctuations proper- 
ties [8]. It has also been shown that the Zeeman split- 
ting leads to the spectrum of the effective single-particle 
Hamiltonian exactly as required by the observed pattern 
of quantization of Hall conductivity p^ . These are in- 
teresting effects associated with the electronic transport 
properties of realistic systems, that we hope to present in 
the future and based on the model we have studied here. 
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